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PREFACE 


■UJU  V  V  'i.1  V’1.’ 


The  models  and  algorithms  in  this  paper  represent  a  synthesis  of 
results  from  different  works  by  several  different  individuals  over  the  past 
few  years.  In  particular,  I  would  like  to  acknowledge  the  contributions 
of  the  following  people:  Dr.  Walter  Chesnut  and  Dr.  Charles  Rino  of  SRI 
International,  Dr.  K.  C.  Yeh  and  his  colleagues  at  the  University  of  Illinois 
at  Urbana,  Dr.  David  Sachs  of  Science  Applications  Incorporated,  Dr.  Roy 
Hendrick,  Mr.  Robert  Bogusch,  Dr.  Dennis  Knepp  of  Mission  Research 
Corporation,  Dr.  Clifford  Prettie  of  Berkeley  Research  Associates,  and 
Dr.  E.  J.  Preoouw  of  Physical  Dynamics  Incorporated. 
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RADIO  WAVE  PROPAGATION  IN  STRUCTURED  IONIZATION 
FOR  SATELLITE  APPLICATIONS  II 


1 .  INTRODUCTION 

This  report  is  an  extension  of  reference  1  which  developed  the 
algorithms  necessary  to  calculate  the  signal  structure  parameters  and  the 
generalized  power  spectrum  which  represent  the  simultaneous  phase  and 
amplitude  effects  of  propagating  electromagnetic  waves  through  media 
characterized  by  a  structured  index  of  refraction.  Portions  of  that 
work  were  restricted  to  structure  spectral  indices  greater  than  or  equal 
to  two.  Also,  the  effects  description  did  not  cover  the  phase  effects 
only  component  of  the  signal.  This  paper  will  reduce  these  restrictions. 

2.  ENVIRONMENT  CHARACTERIZATION 

The  first  step  in  any  propagation  study  is  to  characterize  the 
ionization  or  equivalently  the  index  of  refraction  fluctuations  of  the 
ionospheric  propagation  environment.  Figure  1  illustrates  the  geometry  of 
a  typical  satellite  link.  The  index  of  refraction  fluctuations,  represented 
schematically  by  the  curved  lines  are  typically  elongated  along  the  magnetic 

A 

field,  t  is  a  unit  vector  along  the  magnetic  field  and  is  a  slowly  varying 
function  of  space  since  the  field  lines  are  curved.  The  z  axis  and  the 

A 

z  unit  vector  represent  the  propagation  line  of  sight.  The  transmitter 

A  A 

by  definition,  is  at  z»0  and  the  receiver  is  at  z»z  .  The  r  and  s 

A  ^ 

unit  vectors  complete  with  t  an  orthogonal  coordinate  system  that  is  used 

A  A 

to  define  the  structure.  The  orientation  of  r  and  s  is  chosen  to 
best  represent  any  anisotropy  of  the  index  of  refraction  structure 

A  A  A 

about  the  field  line.  As  with  t,  r  and  s  may  be  slow  functions  of 


3 


position.  In  the  r,  s,  and  t  system  the  index  of  refraction  fluctuations 
are  often  represented  by  a  power  law  power  spectral  density. 


PSD(Kr,Ks,Kt) 


8*3/2An?L  L  L.r  (n)/r(n-3/2) 

a  r  3  i 


(1 


♦L‘K2 


r  r+LsKs+LtKJ)n 
r  r  s  s  t  t 


(1) 


where 


*r*^s’*t  *  spatial  wave  numbers 
L  ,L  ,L.  ■  structure  outer  scales 

T  5  t 
— j 

An^  *  index  of  refraction  variance 
r(n)  *  gamma  function  of  argument  n 


2n  -  2  *  fluctuation  spectral  index 

2-9  10 

Both  in  situ  measurements  and  theoretical  considerations  imply  that 
n<*2  for  ionospheric  structured  environments.  This  spectral  slope  is 
believed  to  be  caused  by  the  existence  of  very  shaip  edged  structures.  The 
fluctuation  power  spectrum  is  considered  a  locally  homogeneous  quantity  as 
discussed  by  Tatavski  11 .  This  means  that  the  parameters  can  be  slowly 
varying  functions  in  any  space  direction  with  respect  to  the  spatial  scale 
in  that  direction. 


The  structure  variation  of  the  index  of  refraction  perpendicular 
to  the  z  axis  dominates  -'e  propagation  effects  while  the  variation  parallel 

A 

to  the  z  axis  enters  only  through  the  strength  of  the  integrated  phase 
variance.  Thus  Equation  1  must  be  transformed  to  a  frame  with  one  axis 

A 

being  the  z  axis.  These  new  axes  are  defined  by: 
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A  A  »  »  A 

X  -  t  X  z/|t  X  t| 


(2a) 


y  ■  x  x  (t  x  z)/|z  x  (t  x  z) J  (2b) 

A 

Two  rotations  suffice  to  accomplish  the  transformation.  First,  the  r  axis 

A  A 

is  rotated  about  the  t  axis  by  angle  $  to  become  parallel  to  the  x 
axis,  *  is  defined  by 

A  A  A  A 

tain*  *  r  x  (t  x  z)  (3a) 

A  A  A 

Cos*  ■  r  •  (t  x  z)  (3b) 

A  A  A 

Next,  the  t  axis  is  rotated  about  the  x  axis  by  -ngle  *  into  the  z 
axis.  *  is  defined  by 

xsin*  »  t  x  z  (4a) 

A  A 

cos*  ■  t  •  z  (4b) 


These  transformations  can  be  simplified  by  defining  new  effective  scale  sizes. 


L2  ■  L“cos2*  ♦  L2sin2*  (5a) 

a  r  3 

L2  -  (l2sin2*  ♦  L2cos2*)cos2*  ♦  L2sin2*  (5b) 

y  r  s  * 

L2  ■  (L2sin2*  +  l2cos2*)sin2*  +  L2cos2*  (5c) 

z  r  s  t 

LXy  »  (L2-l2)cos*cos*sin*  (5d) 

L  ■  (L2-l2)sin*cos*sin*  (5e) 

xz  r  s 

L  »  (L2-L2sin2*-L  cos2*)cos*sin*  (5f) 

yz  tr  s 
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The  final  result  for  the  local 1/  homogeneous  index  of  refraction  power 
spectral  density  is 


8ir3/2An?L  L  L  r(n)/r(n-3/2) 

POL.K.10  *  — vryi"  - - n 

y  (l+L;iT-*i/ir+Lzir+2L  K  K  +2L  K  K  +2L  K  K  )n 

x  x  y  y  z  z  xyxy  xzxz  y:  y  r 


DERIVATION  OF  SIGNAL  STRUCTURE  PARAMETERS 


The  required  signal  parameters  are  derived  from  the  differential 

phase  spectrum  which  is  Equation  6  with  K  set  to  zero  and  multiplied  by 

2  1 
K  ,  wheie  K  is  the  carrier  frequency  wave  nuuber.  The  differential  phase 

spectrum  is 


dP^(Kx>Ky)  8wV2IC2An2Lrl5Ltr  (n)  /r  (n-3/2) 

dI  "  O'^’VxV" 


The  mean  square  integrated  phase  fluctuation  is 


’♦ '  Jil  “a? 


where 


d%  _  iAoi-ii  .2“i. 
IT  '  'rfa-Vri  KtaiL, 
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Cio) 


j^prsronr 


\N 

v« 


w  %r‘.  ’V.V 


■■-  •  • « .  %  *. 


IV.*.'  *.  .V 


k*,1 

t\  ( 


L  L  L 
,  _  rLsbt 

LI  ,757 i  ; r  J, 


(L  L  -L  ) 
x  y  xy 


I*z  is  redefined  as  the  effective  z  axis  scale.  Finally 


<jy*,.y  _  *£  *■  (L»Ly-Lw)h  tn~u 


(1+L  K %L  K %2L  K  K  ) 
v  x  x  y  y  xy  x  y' 


(11) 


Equation  9  requires  n  greater  than  1.5  which  is  an  artifact  of  the  spectral 
form  chosen  for  the  index  of  refraction  fluctuations.  Equation  11  requires 
n  greater  than  one  if  doZ/dz  is  suitably  redefined  consistent  with 
Equation  8.  The  following  developaent  applies  to  1.2sn*4  where  a  well 
defined  do^/dz  is  assumed. 


Equation  11  can  be  Fourier  transformed  to 


d^(x.y)  da* 
dz  *  dz 


n-l„ 
P  K. 


n-l<P) 

2n'Zr(n-l) 


(12a) 


{'-V° 


... 


(12b) 


where 


2  '.y^L^y 

2  2  2 
L*L  -l 
x  y  xy 
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m  ■  2,  n>2 


m  =  2n-2,  1.2sns2 

Bn  =  constant  of  order  unity  (see  Appendix  B) 

Equations  8,  11,  and  12  are  the  fundamental  equations  from  which  all  of 
the  required  propagation  quantities  can  be  derived. 

The  mean  square  log  amplitude  fluctuation  is  developed  in 
Appendix  A  where  the  range  of  application  has  been  extended  from  that  in 
Reference  1. 


~2 

nr  “ 


CZ t  .  2 

I  fdz  (n“ 


6  2  2  4 

8a°  8(a;-cpc^ 

1)£(n)  ^  3M(z)  +  M(zl 


where 


2^4^-! 


+  lo™  + 


16ciC2 


(z)a: 


2-2n 


2n-2 


6-2n 


lec^V'^zlCc^211^'211)  ,  n^3 


(13) 


lec^V'^zDlnW,  n*3 


M(z) 


(L^L2)(zt-z)z 

Kjt(Lxly-Lx^ 


k 

aj  *  rain  (c^CjM^fz)) 
v 

a2  *  max  (Cj.^M  (z)} 

Cj  =  0.5 
c2  =  0.841 

f(n) .  ™ ' '1 
f'  (n)  *  1.1  -  max  [0.0,0.5(n-2 .  4)] 


The  perpendicular  signal  decorrelation  as  developed  in  Reference 
1  already  applies  over  the  required  range  of  spectral  indices. 


where 


(14) 


Vx-x.V 


Vy  *  y  *  v 
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V  V  MM\'  It'-  <’l>  '  »  ’  *  ’  « 


V.  *  transmitter  velocity 
tr 

V  =  index  of  refraction  structure  velocity 
Vre  =  receiver  velocity 

B(n,cr2)  -  minimum  [^1,  (-0.34n2+2.51n-2.00)  (o^Bn)1/mJ 

The  remaining  signal  structure  parameters  are  derivable  from  the 
mutual  coherence  function  of  U(z,r,K)  where  the  propagating  carrier 
signal  is  expressed  as 


-|2 


E(z 


,r,K)  . 


(15) 


The  mutual  coherence  function  for  U  is 


where 


The  equation  for  G  is 


dG 

Hz 


.  n)  -  U^z.rj 

JCjM* 

»r2»K  2^ 

-k2j/k2 

V71 

Ji2G 

vd2c\ 

x  dG  y 

M-") 

v) 

’  z  dx  ’  z 

\iL 

1 

y°l 

<Vx,y) 

K 

d-n) 

2  I  dz 

dz  1 

11 

(16) 


(17) 


I  "  •  *  4  ^  • 


»?.  -.■  r."  v" •>  -  '  -  *.*"*.•  V  V  "  • 


First,  an  approximate  parallel  signal  decorrelation  time  and 
angle  of  arrival  distribution  will  be  developed.  Let  n  be  zero.  Also, 

A  A  A 

let  us  define  a  set  of  fixed  axis,  u  and  v,  perpendicular  to  the  z  axis 

AAA  A 

such  that  u  x  v  ■  z.  At  each  point  on  that  z  axis,  the  rotation  angle 

A  A 

between  the  x  and  the  u  axis  is  defined  by 

A  A  A 

z  sin  0 ■  x  x  u 

A  A 

cos  0  *  X  •  u 


L2  ■  L2cos20+L2sin20-2L  sin8cos0 
u  y  x  xy 


L2  *  L2cos20+L2sin20+2L  sin0cos0 
v  a  y  xy 


Luv  ■  (Ly-L2)sin0cos6-Lxy(sin2e-cos20) 


The  formal  solution  to  Equation  17  for  o^>l  is 


G(zt,u,v) 


/>  iai .  L 


The  exponential  is  not  factorable  so  we  approximate  Equation  18  by 


G(zt,u,v)  »  e 


-(C  u2+C  v2-2C  uv)"^2 
v  u  V  uv  * 
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where 


f  4  4)’ 


>/2 

dz 


i2/m 


.  2 

2 

m/2 

ft  <% 

1  ^  " 

A  V  v  4 

dz 

VV  LxLy"L^y 

2/m 


(20a) 


(20b) 


CUV  ■  si«"«V  l^.vl2/" 


(20c) 


SW 


■p  *4 


2  l«tol 
tf  lV-L? 


W 


x  y  xy 


m/2 

sign(Luv)dz 


(20d) 


Equations  19  and  20  are  approximate  generalizations  of  Equations  25  and 
26  of  Reference  1.  Equations  19  and  20  agree  with  Equation  18  when  n 
is  greater  than  or  equal  to  two  or  if  Lu,  Ly,  and  Luv  are  proportional  to 
each  other  along  the  z  axis.  Thus,  Equations  19  and  20  are  always  reason¬ 
able  approximations  for  satellite  applications  independent  of  n. 


Next,  let  us  rotate  the  u-v  coordinates  to  remove  the  cross  term 
in  G.  Let 


2C 

tan  (2e)  - 

u  v 
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p  ■  u  cos  c  -  v  sin  c 


Then 


q  ■  «  sin  c  ♦  v  cos  c 
cp  -  7 
c,  *  7 


RUt.p»q) 


-eyV^2 


(21) 


The  parallel  signal  decorrelation  time  is  defined  as 


T*  (Cp/3+C^3)3/2|(Vre-V5t)*z|  (22) 

The  final  decorrelation  time  is 

tq  -  minimum  (t^.t,)  (23) 

The  angle  of  arrival  distribution  is  approximated  by 


WV 


-tcpp2*cq<i2)*/2  »p^> 

0  0 


ie  Kq 

e  M 

(24) 
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where 


VVV  ’ 1 


(25) 


When  ■  equals  two.  Equation  24  is  equal  to  Equation  36  in  Reference  1. 
The  measures  of  the  angular  scatter  are  defined  by 


IT 

(26a) 

% 

IT 

(26b) 

For  m  equal  to  two.  Equation  26  is  the  aean  square  scattering  angle.  For 

a  less  than  two,  this  interpretation  is  not  strictly  valid.  Nevertheless, 
2  2 

<rfi  and  o.  rent  in  good  aeasures  of  angular  scatter. 

P  q 

The  generalized  power  spectrins  can  be  derived  by  solving 
Equation  17  using  Equation  12  and  the  delta  layer  approximation. 

d04  f  t<io* 

JT -  dCz-Zj)^  -gj  dz  (27) 


To  account  for  time  decorrelation,  the  argument  for  Equation  12b  is 
Modified  by  incorporating  the  effective  velocity,  Vx,  and  a  time 
displacement  argument 


(28) 
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The  cross  ten  has  been  eliminated  in  anticipation  of  the  two  final  cases: 
one  dimensional  structure  perpendicular  to  the  line  of  sight  and  isotropic 
structure  around  the  line  of  sight.  The  effective  velocity  is  assumed 
parallel  to  the  %  axis.  After  some  algebra  we  get 


16 


where 


2 

o^R  ■  Rayleigh  phase  variance  (see  Appendix  A) 

f  ■  scintillation  frequency  spread 
t  ■  delay 


. ..  hw^r 

— j 


9x**ym  energy  arrival  angle  coordinates 


The  final  generalised  power  spectxrm  it  ca.  nJLated  by  integrating  the 
product  of  Equation  29  and  the  transmitter  and  receiver  antenna  gain 
functions  over  the  energy  angle  of  arrival  coordinates. 


rn(f.r) 


x,  v*.*>  vv  vG’Nr  vsr  v 


(30) 


where  is  the  root  aean  square  energy  angle  of  arrival  with  the 
transmitter  and  receiver  interchanged -  Calculations  shows  that  the 
functional  fora  in  Equation  30  n*?d yell  for  thick  layers  as  long  as 
r  ,  f  ,  ofi,  and  o£  are  choser.  properly .  Appendix  D  derives  a  slightly  nore 
general  fora  of  Equation  30  and  *  useft..*:  approximation  for  Equation  30  for 

spectral  indices  greater  than  or  equal  to  two.  The  Rayleigh  phase  variance, 

2 

'>R*  was  introduced  in  Reference  1  to  distinguish  between  that  portion  of 
the  phase  variance  that  accounts  for  the  Rayleigh  component  of  the  signal 
and  that  portion  that  represents  a  phase  effects  only  randoa  signal 
component .  f£  and  f£  as  defined  above  apply  to  the  delta  layer.  A 
reasonable  extension  for  thick  layers  is  available  by  applying  a  scale 
factor,  H^,  to  Equation  38  in  Reference  1. 
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where 


fj*  ’  7?1  I  7*1?  I 


f1 '  do2  /2i2L  2  \ 

V,,)  mr'^bfir) 

<  \  X  y  xy/ 

wt,,)  'b 

» il  “»^J 


and  t!  is  redefined  as 


Equations  29  through  32  provide  a  reasonable  representation  of  frequency 
selective  effects.  For  n  greater  than  or  equal  to  two,  the  results  are 
quite  accurate.  For  smaller  n,  the  results  are  less  accurate.  The  n 
equal  to  3/2  case  was  examined  in  detail.  The  inaccuracy  in  Equations 
29  through  32  were  quite  acceptable  including  the  case  where  both  the 
transmitter  and  the  receiver  were  in  the  turbulent  medium.  The  errors 
decrease  as  the  layer  thickness  becomes  comparable  or  saaller  than  the 
free  space  propagation  distances. 

The  frequency  selective  bandwidth  is  still  defined  by: 


4. 


THE  PHASE  EFFECTS  ONLY  SIGNAL  COMPONENT 


(33) 


In  the  previous  section,  it  was  mentioned  that  a  portion  of  the 
integrated  phase  spectrum  determines  a  phase  effects  only  signal  component. 
This  signal  component  can  be  represented  as 


(34) 


where  4^(t)  represents  possible  large  scale  trends  in  the  phase  and  g(f) 
is  a  zero  mean  normally  distributed  variable  representing  the  random  phase 
component.  The  total  signal  is  the  product  of  Equation  34  and  the 
Rayleigh  component.  The  variance  of  g(f)  is 


PTTOFT  -  «  (f-f')P^(f) 


(35) 
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where  g(>£)  equals  g*(f).  An  approximate  expression  for  the  phase  power 
spectrin,  P^(f) ,  can  be  derived  using  the  delta  layer  approximation. 

Let 


*  *  f  v 

t 

(36a) 

y  *  Ir  V 

(36b) 

From  Equation  12  and  the  delta  layer  approximation,  the  phase  time 
correlation  function  is 


R.U)  n 

*  2n”*r(n-l) 


$  °A  v) 
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Also  let 


[<•>/'■  isf 
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The  phase  frequency  spectrin  is 
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The  phase  power  spectrum  should  be  truncated  at  fR  where 


m 

2  /p4(f)df  -  max(oJR  0.025) 


(40) 


The  frequencies  greater  than  fR  should  be  ignored  because  either  they 
contribute  to  the  Rayleigh  signal  component  or  they  contribute  little  to 
the  total  integrated  phase  variance.  Prom  Equations  39  and  40 

1 

1ST? 

(4!) 


,53'i- 

2irr  ir^(n- 

on  l 


r  (n-%) 


1)  r(n-l)max(0 


♦R 


2  2 

In  Equation  41,  fR  is  set  to  zero  when  o^R  equcus  o^.  Equation  39  will 
be  equal  to  or  greater  than  the  actual  phase  spectral  power.  The  actual 
phase  spectrum  for  a  given  n  will  be  shallower  in  slope  than  above 
because  of  geometric  or  physical  non-global  homogeneities  in  the  index 
of  refraction  statistics.  The  two  will  agree  at  approximately  f  oqual  to 
fR,  thus  the  actual  spectrum  will  fall  below  Equation  39  for  frequencies 
less  than  fR. 


In  applications  where  the  structured  index  of  refraction  is  due 
to  ionization  structure.  Ihe  phase  effects  only  phase  spectrum  is 
usually  represented  as  a  total  electron  content  power  spectrum. 

PTECtf)  ■ 

where  r  is  the  classical  electron  radius. 

0 


(42) 


The  development  in  this  and  the  previous  section  assumed 
that  n  was  a  globally  homogeneous  quantity.  Proposed  algorithms  for 
lifting  this  restriction  are  presented  in  Appendix  C. 

5.  SUMMARY. 

The  proceeding  sections  extended  the  results  of  Reference  1  to 
permit  the  prediction  of  the  effects  of  structured  index  of  refraction 
on  electromagnetic  waves  to  spectral  indices  of  less  than  two.  Ihe 
phase  effects  only  phase  spectrum  was  developed,  which  along  with  the 
Rayleigh  component  constitute  a  complete  statistical  description  of  a 
scintillated  signal.  This  work  and  that  in  Reference  1  form  the  basis 
of  a  signal  structure  specification  for  satellite  C3  applications 
proposed  in  Reference  12. 
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APPENDIX  A 

CALCULATION  OF  THE  MEAN  SQUARE  LOG  AMPLITUDE  FLUCTUATION 
AND  THE  RAYLEIGH  PHASE  VARIANCE 


The  equation  for  the  aean  square  log  amplitude  fluctuation  is 


7 


(A-l) 


where 


dp»t>V  .  4,fL>lr-L^1/2tn-1) 

dz  (i+lV*it+2l  n)n 

k  x  x  y  y  xy  x  yJ 


(A- 2) 


Equation  A-l  using  A-2  is  not  integrable  but  it  can  be  reasonably- 
approximated.  Let 


Kx  ■  ic  cos  6 


Ky  ■  K  sin  0 
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2  2  4 

sin  (a  j  -  a\  a<c2 


-  c2,  a>c2 


The  differential  phase  spectrum,  in  tens  of  a,  is  assumed  constant  for 


acCjM^z) 


and  power  law  for  larger  values.  Pinal ly  let 


ax  -  minCcj.CjM^Cz)) 


a2  "  Biax 


(c^Cjl^Cz)) 


The  final  approximation  for  1 .2<ns4  is 
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This  range  assumes  that  %  is  well  defined  for  n  less  than  or  equal  to 
1  5.  dz 

With  the  following  choices  for  c^,  c^,  and  f(n),  the  erroi  in  Equation 
A-3  is  less  than  thirty  percent. 


Cj  «  0.5 
c2  ■  0.84 


f(n) 

f’(n) 

The  correction  factors  apply  for  all  M(z)  and  n  between  1.2  and  4.  Tbey 
also  apply  for  any  value  of  L^Ly . 

Equation  A-3  is  consistent  and  a  bit  nore  accurate  than  the  x2 
dovelopaent  in  Reference  1.  That  result  is  still  often  preferable  over  its 
range  of  application  because  of  its  simplicity. 


The  phase  variance,  for  snail  M(z),  is  divisible  into  the  Rayleigh 
phase  variance  and  the  phase  effects  only  portion.  The  Rayleigh  phase 
variance  is  defined  as 
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'♦R  ^ - J^T T 


ac  is  chosen  using  Equation  A-3  to  determine  that  portion  of  the  integral 
~7  “T 

that  contributes  xc  to  x  .  This  condition  is  approximated  by 
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Xc  is  chosen  as  0.1.  The  procedure  to  calculate  o.f  is  relatively  straight 

-5  2  ~T~ 

forward.  First,  if  x  <0.1,  then  a  *»0.  End  procedure.  If  lat  ,  then 
22  c' 

o^R  ■  0^.  End  procedure.  Otherwise  calculate 


r*t  2 
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Next  calculate 

1 j  -  (I2c2~2n  ~  IicJ_2n)/(6-2n),  m*3 
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It  lj<0  and  o^/2<Xc,  then  <j^r«0.  End  procedure.  If  I3<0  and  o^/2>Xc> 

2  2 

then  o^R»c^.  End  procedure.  Otherwise  calculate 
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If  X  -  X  <  0,  then  calculate 
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By  twice  differentiating  Equations  B-l  and  B-2,  it  is  easy  to  show  for  any 
p  less  then  about  0.1 


Im-ISS. l~  k2  co»  (.k)  dk 
*  v’r(n-l)  l  (l*k2)n-’ 

Equation  B-3  can  be  approximated  by 


b  s,  -Maa . 

n  cnir*r(n-l) 
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where  f(n)/(n-l)  is  a  correction  factor  to  correct  for  the  approximation 
of  the  integral  in  Equation  B-3.  Now 
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In  the  limit  as  e  becomes  small 
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By  using  Equation  B-6  and  requiring  continuity  of  Bn  at  n  equal  to  two, 
the  correction  factor  can  be  evaluated 

f(n)*  (n-1)'0*5,  n>2  (B-7) 

For  n  less  than  two,  let  us  differentiate  Equations  B-l  and  B-3 
once .  Then 


cnir^(n-l)r(n-l) 


x  sin[(l+e2/p2)x1  dx 


,  2  2.n->j 
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where  x  equals  kp.  Equation  B-8  can  be  approximated  by 
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sin(x)dt 
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where  it  is  assumed  that  p  is  greater  than  or  equal  to  Ge.  Now 
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As  expected,  Equation  B-10  is  relatively  insensitive  to  p  which  can  be 
thus  thought  of  as  a  fit  parameter.  By  choosing  p  equal  to  e  and  fitting 
f(n).  Equation  B-10  can  be  made  continuous  with  Equation  B-5.  For 
small  e,  it  can  be  shown  that 
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where 
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g(n)  -  j  ,  n»y 


g(n)  -  r(3-2n)  sin  ^-3-^?^,  l<n«| 


Using  Equetion  B-ll  and  the  continuity  requirement  at  n  equal  to  two 
find 


f(n)  -  (n-1) 
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The  final  approximation  for  Bn  is 
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where 
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Equation  B-13  is  accurate  to  within  twenty  percent  for  the  specified 
variable  ranges. 

In  principle,  the  inner  to  outer  scale  size  ratio  is  not  an 
isotropic  quantity.  Fortunately,  Equation  13  is  insensitive  to  the 
precise  value  of  e,  hence  any  reasonable  estimate  will  do.  To  avoid 
underestiaating  effects,  the  aaxiaua  reasonable  value  of  e  should  be  used. 
If  this  procedure  is  inadequate,  then  the  anisotropy  can  be  handled  by 
defining  Bn  as  the  non-c  dependent  portion  of  Equation  B-13  and  redefining 
the  outer  scales.  For  example,  if  n  equals  two,  then  the  following 
definitions  can  be  used. 
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1  ,  and  lt  are  the  inner  scales.  These  new  quantities  substitute 
,  L^,  and  BR  in  all  of  the  equations  in  this  report. 
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APPENDIX  C 


LOCALLY  HOMOGENEOUS  SPECTRAL  INDICES 


The  developments  in  Reference  1  and  the  present  paper  have 
usually  assumed  that  the  fluctuation  spectral  index,  or  equivalently,  n, 
is  a  globally  homogeneous  variable.  This  appendix  will  propose  procedures 
to  lift  this  limitation.  It  is  assumed  that  the  reader  is  familiar  with 
the  other  portions  of  this  paper  including  the  preceding  appendices. 

Let  us  start  with  the  mean  square  log  amplitude  fluctuation  and 
the  Rayleigh  phase  variance.  Equation  A-3  already  applies  for  variable  n. 
In  preparation  for  calculation  of  the  Rayleigh  phase  variance.  Equation  A-3 
must  be  divided  up  into  separate  contributions  from  discrete  values  of  n. 
Let 
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If  <0,  then  End  procedure.  Otherwise  solve  for  ac  in  the 

following  equation. 
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The  next  parameter  is  the  perpendicular  signal  decorrelation  tine, 
t.  .  First,  let 
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ta  then  satisfies  the  following  equation. 
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We  also  define  the  effective  time  index 
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Now  calculate  au  and  ay  where 
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The  angular  index  is  defined  as 
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The  above  coefficients  are  then  used  in  Equation  19  with  a0  substituting 
for  a.  Thus 
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Equations  21  through  26  then  follow.  Equation  30  remains  a  reasonable 
choice  for  the  generalized  power  spectnna  where  the  index,  a,  is  redefined 

as 


a  -  ainiaua  (a0,  aT)  (C-15) 

where 

ax  -  ainiaua  (2,  2nx-2)  (C-16) 

Equation  C-15  imposes  a  worst  case  functional  fora  for  the  generalized 
power  spectrum  with  regards  to  both  the  frequency  and  delay.  For 
situations  where  frequency  selectivity  is  not  significant  it  is 
best  to  use  aT  for  a  to  avoid  ov eras tiaa ting  the  fade  rate  effects. 
Conversely,  if  frequency  selective  effects  doainate,  aQ  is  probably  the 
best  choice.  If  both  effects  are  operative,  then  Equation  C-15  is  advised 
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so  ss  to  not  underestimate  the  effects.  In  any  event,  non-globally 
homogeneous  spectml  indices  imply  soee  compromise  in  the  fidelity  of 
the  generalized  power  spectrum  form  so  some  conservatism  is  good. 

Equation  31  through  33  remain  reasonable  with  me  replacing  m 
in  Equation  31.  Similarly  Equations  39  and  41  can  be  used  with  n^  and  mT 
replacing  n  and  mfi  respectively. 

The  choices  for  n^  and  An^  depend  on  the  range  of  n  and  the 
required  resolution.  If  n  varies  continuously,  then  several  values  of  n 
are  necessary  with  the  finest  resolution  for  small  values  of  n.  Another 
possibility  is  where  n  has  a  number  of  discrete  values  corresponding  to 
different  physical  processes  in  different  regions.  In  this  case,  An^ 
is  chosen  as  some  arbitrarily  small  value.  The  formalism  does  not  allow 
equal  to  three,  but  any  choice  arbitrarily  close  to  three  is  adequate. 

This  completes  the  formalism  for  locally  homogeneous  spectral 
indices.  Unfortunately,  this  formalism  is  significantly  more  complex  and 
costly  to  implement  over  that  for  a  globally  spectral  index.  This 
formalism  also  increases  the  degree  of  approximation  particularly  with 
regards  to  Equation  C-14,  the  functional  form  of  the  generalized  power 
spectram,  f£,  and  the  functional  fora  of  the  phase  effects  only  phase 
power  spectrum.  These  results  still  provide,  however,  reasonable 
estimates  of  the  signal  structures  and  their  resulting  effects. 
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APPENDIX  D 

THE  TWO  FREQUENCY  MUTUAL  CORRELATION  FUNCTION 
AND  ANTENNA  EFFECTS 


To  derive  the  various  results  of  this  appendix,  the  most 
general  equation  for  the  two  frequency  autual  correlation  function  is 
required.  Let  that  function  be 


G(i,x,y,X,Y,n)  »  U(z,r2,K2) 


where  the  orignial  propagating  signal  is  expressed  as 
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The  equation  for  G  is 
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Equation  17  results  from  Equation  D-3  by  neglecting  variation 
of  R^(x,y,X,Y)  with  respect  to  X  and  Y  and  by  assuming  no  initial 
dependence  of  G  on  X  and  Y. 


Equation  D-3  can  be  used  along  with  the  delta  layer  approximation 
and  considerations  from  antenna  theory  to  get  the  generalized  power  spectrum 
including  antenna  effects.  Let  us  assume  that  the  receiving  and  trans¬ 
mitting  antennas  are  circular  aperatures  with  gains  of  G  and  G', 
respectively.  In  tens  of  the  antenna  aperature  diameters 


where  X  is  the  propagation  wavelength. 


(D-4a) 


(D-4b) 
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Assuming  isotropic  fluctuations  about  the  propagation  line  of  sight 
wd  the  generalized  power  spectrum  is 


where 


s2  -  G  |^Tofa0+ex)2  ♦  p\Roe+0y)2j 

s'2  ■  G*  £(2Vrofo£+e^) 2  +  ^irRo^+e^J 
0x,0v  *  **®ceiving  antenna  pointing  error 

®y  *  transmitting  antenna  pointing  error 
Jj(s)  ■  Bessel  function 
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Rayleigh  phase  variance 


*2 

<7g  *  energy  angle  of  arrival  variance  at  transmitter 

m  ■  minimum  of  two  and  spectral  index 

For  one  dimensional  structure  transverse  to  the  propagation  line  of 
sight,  insert  5  (M)  6  (R)  at  the  end  of  Equation  D-S  and  integrate  over  R 
and  M.  The  gain  factors  and  the  free  space  propagation  loss  have  been 
left  in  Equation  E-S.  Usually  these  terms  are  handled  separately. 


By  dropping  these  terms,  assuming  a  spectral  index  greater  than 
or  equal  to  two,  zeroing  the  pointing  errors,  and  approximating  the  Bessel 
functions  terms  by  a  gaussian  function,  the  following  generalized  power 
spectrum  can  be  calculated. 
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and  the  Bessel  function  terms  were  approximated  by 
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The  functional  form  of  the  power  spectrums  is  identical  to  that  in  Reference 
1.  Thus,  the  results  in  Reference  1  and  12  can  be  used  after  scaling  tq 
and  fQ  and  accounting  for  the  power  loss.  The  scaling  of  tq  and  f  is 
generally  beneficial  to  systems.  Unfortunately,  the  power  loss  generally 
negates  much  of  the  gain  unless  there  is  a  large  margin. 
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